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Boltzmann-Gibbs statistics fails to study the systems having the condi- 
tions (i) the spatial range of the microscopic interactions are long-ranged, (ii) 
the time range of the microscopic memory is long-ranged and (iii) the system 
evolves in a (multi) fractal space-time. These kind of systems are said to be 
nonextensive and a nonextensive formalism of statistics must be needed for 
them. Recently a generalized thermostatistics is proposed by C. Tsallis to 
handle the nonextensive systems and up to now, not only the generalization 
of some of the conventional concepts have been investigated but the formalism 
has also been properous in some of the physical applications. In this study, 
our effort is to introduce Tsallis thermostatistics in some details and to give 
a brief review of the magnetic systems which have been studied in the frame 
of this formalism. 

I. INTRODUCTION 

The aim of Statistical Thermodynamics is to establish a bridge between microscopic and 
macroscopic behaviours of all types of systems in Nature. Standard Boltzmann-Gibbs (BG) 
statistics (or in other words extensive statistical thermodynamics) is essentially derived from 
the standard Shannon entropy 

w 

Si = -k B J2pi\ogpi (1) 

i=l 
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which is an extensive and concave quantity (the sub index 1 of S will be defined afterwards). 
Although BG statistics provides a suitable tool which enables us to handle a large number 
of physical systems satisfactorily, it has two basic restrictions : 

• the range of the microscopic interactions must be small compared to the linear size of 
the macroscopic systems (short-range interactions) 

• the time range of the microscopic memory must be small compared to the observation 
time (Marcovian processes). 

In the case of a breakdown in one and/or the other of these restrictions, BG statistics fails. 
Here, "to fail" is used to imply the divergences of the sums or integrals in the expressions 
of the relevant thermodynamical quantities. 

These kind of violations are met for a long time in gravitation [1], magnetic systems [2], 
anomalous diffusion [3] and surface tension problems [4]. On the other hand, recently C. 
Tsallis have pointed out that [5], although not yet clearly identified, the same or analogous 
type of problems might be present in long-range Casimir-like systems [6] , in granular matter 
[7] and two-dimensional turbulence [8,9]. More precisely, the situation could be classified in 
a general manner as follows [5]: 

• For an Euclidean-like space-time, if the forces and/or the memory are long-ranged, as 
far as we are interested in an equilibrium state, the BG statistics is weakly violated, 
therefore BG formalism can be used. On the other hand, whenever a meta-equilibrium 
state is considered, the BG statistics is strongly violated, hence another formalism must 
be needed. 

• For a (mult i) fractal space, BG formalism is strongly violated again and a new formal- 
ism is needed. 

The way out from these problems seems to be Nonextensive Statistical Thermodynamics 
which must be a generalization of the BG statistics in a manner that allows a correct 
description of the nonextensive physical systems as well. 
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Nonextensive formalisms are much in vogue nowadays in Physics and keep growing 
along two apparently different lines, namely Generalized Tsallis Thermostatistics (GTT) 
and Quantum Group-like Approaches (QGA). Although they seem to be very different from 
each other, in a recently published paper [10], it is noticed that there exist a connection 
between them. Since this is not the scope of the present paper, we're contented to imply 
that there have been attempts of finding some verifications of this connection or establishing 
new relations by investigating the physical systems within both GTT and QGA [11,12]. 

From now on, we focus our attention to GTT and the following Section is devoted to a 
review of the formalism. 



II. GENERALIZED TSALLIS THERMOSTATISTICS 

GTT, which has been proposed in 1988 [13] as a possible theoretical tool for discussing 
nonextensive physical systems, basically relies upon two postulates: 

1) The generalized entropy is assumed to be defined as 

1 _ V w n q 

S q = k _i ( 2 ) 

where k is a positive constant, q G Sft, W is the total number of microscopic configurations and 
{pi} are the probabilities of the microstates of the systems. This expression is a nonextensive 
quantity and transforms to the standard extensive Shannon entropy given by Eq.(l) if and 
only if q — 1. 

2) The g-expection value of an observable O, whose value in state % is O iy is given by 

w 

<0> q =Y,PiO l (3) 

i 

It is clear that the validity of these postulates is confirmed (or rejected) by comparisons to 
the experiments and by the conclusions which they yield. It is easy to verify that GTT has 
the following properties: 

• S > for any arbitrary set {p^ and for any value of q (positivity). 



• For microcanonical ensemble (i.e. pi = 1/W , Vi), S q attains its extremal value: 
S q = (W l ~ q - 1)/(1 - q) (equiprobability). 

• S q is concave for q > and convex for q < (concavity). 

• The optimization of S q under the constraints 

w w 

J2Pi = l and U q = Y,Pi £ i, ( 4 ) 

i i 

leads to the canonical probability distribution 

^ = ^[l-(l-^] 1/(1 ^ (5) 

where the generalized partition function is defined to be 

w 

Z^^tl-a-g)^] 1 ^ (6) 

i 

and = 1/kT. 

• If A and B are two independent systems {i.e. pf^ B = pfpf), S q obeys the following 
additivity rule: (pseudo-additivity) 

S q (A UB)= S q (A) + S q (B) + (1 - q)S q (A)S q (B) (7) 

which implies S q is superadditive (entropy of the whole system is greater than the sum 
of its parts) for q < 1 and subadditive (entropy of the whole system is smaller than 
the sum of its parts) for q > 1. 

• Under general conditions [14], dS q /dt is non-negative for q > 0, vanishes for q = 1 and 
non-positive for q < 1. (//-theorem). 

• If we partition W microstates into two subsets x and y having respectively W x and 
W y microstates (W x + W y = W), corresponding probabilities can be written as 
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Pi,P2, ...,pw x +i an d Pw x +ii ■■■iPw- Thus, it is straightforward to verify the following 
additivity rule: (Shannon additivity) 

S q ( Pl , ..., PW ) = S q ( Px , Py ) + plS q ( * P -^) + plS q f^±i, P A . (8) 

\Px Px J \ Py Py J 

• In this formalism, the Legendre-transform structure of Thermodynamics is invariant 
for all values of q [15]. 

The last property is very important and deserves a further comment. As it has been pointed 
out in [16], this property indicates that the entire formalism of Thermodynamics can be 
extended to be nonextensive without loosing its Legendre-transform structure. In addition 
to this, amongst a large number of entropic forms, Tsallis formalism seems to be the unique 
which preserves the standard Legendre-transform structure. 

From the year 1988 up to present days, numerous concepts of statistical thermodynamics 
have been accomplished to generalize in the frame of Tsallis formalism. Amongst them, the 
specific heat of the harmonic oscillator [17], one-dimensional Ising model [18,19], the Ehren- 
fest theorem [20], the von Neumann equation [21], quantum statistics [22], the fluctuation- 
dissipation theorem [23], Langevin and Focker-Planck equations [24], the Bogolyubov in- 
equality [25], classical equipartition theorem [26], paramagnetic systems [27], infinite-range 
spin-1/2 Ising ferromagnet [28], Callen identity [29], mean-field Ising model [30,31], Planck 
radiation law [32,33], quantum uncertainty [34], anisotropic rigid rotator [35], Haldane ex- 
clusion statistics [36], Bose-Einstein condensation [37], the generalized transmissivity for 
spin-1/2 Ising ferromagnet [38], self-dual planar Ising ferromagnet [39] and localized-spins 
ideal paramagnet [40] could be enumerated. 

GTT has also been successfully used to overcome the failure of BG statistics in some of 
the physical applications. The establishment of finite mass for the astrophysical systems in 
the frame of polytropic structures [41], the calculation of the specific heat of the unionized 
hydrogen atom [42] , the derivation of Levy-like anomalous diffusion [43-45] , the construction 
of a comprehensive thermodynamic description of d = 2 Euler turbulence [9] and solar 
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neutrino problem [46] are the examples that could be mentioned herewith. A detailed 
review of the subject can be found in [47], and the investigation of the formalism from the 
mathematical point of view is now available in [48]. 

After introducing GTT in some details, the remaining part of this paper includes a brief 
review of the magnetic systems which have been studied within this formalism. 

III. REVIEW OF MAGNETIC SYSTEMS IN THE FRAME OF GTT 

Up to now, some of the magnetic systems have been investigated within GTT and in this 
Section our goal is to review these systems briefly rather than to discuss them in details. 

A. Ising Chain 

The investigation of the Ising chain has been the first example of magnetic systems 
discussed within GTT. In his first paper [18] related to the subject, Andrade has evaluated 
the partition function, the internal energy and the specific heat of the Ising chain by making 
use of the probability distribution given in [13]. He has also discussed the extensivity of the 
thermodynamical quantities and the influence of the ground state energy. In the calculations 
a parameter r, defined asr = l/g — 1 with r = 1, 2, ... , has been introduced. Since q — > 1 
if r — > oo, the increasing values of r provide a comparison to the standard values of BG 
statistics. Such kind of comparison has been attempted recently [49]. Although the internal 
energy and the specific heat could be obtained only for the values of r up to 3 because of the 
cumbersome nature of the expressions, it has been sufficient to make a comparison to the 
conventional expressions of the Ising chain. Other two interesting results of Andrade's work 
are those (i) for the negative energy levels, the expressions of the specific heat deviate from 
the expected behaviour for intensive quantities in BG statistics, namely they are proportional 
to 1/N 2 in the limit of large number of spins and (ii) for the positive energy levels, the 
specific heat presents oscillations. In his second paper on the same subject [19], Andrade has 
reanalyzed the behaviour of the Ising chain within GTT by making use of the probability 
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distribution given in Eq.(4), which is now known as the correct expression (for a discussion 
of this situation see [15]). 

B. Infinite-range Ising Ferromagnet 

In order to define infinite-range Ising ferromagnet, let us write down the Hamiltonian 

H = -2 £ J ijSi Sj ( Si = =Fl) (9) 

where 

Jij = 4fs ( J>0 ; d + s>o) , (io) 

r ij 

Tij is the distance between sites i and j and the sum runs over all distinct pairs of sites on 
a d- dimensional simple hypercubic lattice. Infinite-range Ising ferromagnet can be obtained 
in the case of d + 5 — 0. For the internal energy of this system in the frame of BG statistics, 
it can be written at T = (T being the temperature) 

At long distances this sum can be replaced by an integral which diverges for 5 < 0. This 
fact has first pointed out by Hiley and Joyce [2], and it becomes obvious that BG statistics 
fails for long-range Ising ferromagnet when 5 < 0, hence a nonextensive formalism, such as 
GTT, must be needed. 

Nobre and Tsallis have numerically investigated the infinite-range Ising ferromagnet 
within GTT [28]. They have evaluated the specific heat of the system numerically for 
q — 1, q < 1 and q > 1. Furthermore, they have exhibited that the thermodynamical limit 
is well defined. The most remarkable two conclusions of this paper can be quoted by the 
words of Nobre and Tsallis : 

• In complete analogy with the well known Ci/Nks vs. ksT/N curves within BG 
statistics (q = 1), the curves N 2 C q /2^-^ N k vs. kT/N 2 J within q ^ 1 statistics tend 
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to a numerically well defined thermodynamical limit; this is the first time that the 
existence of such a limit is exhibited for an interacting model 

• Analogously with the Landau-like phase transition which is known to exist for q = 1, 
a nontrivial divergence, at a finite rescaled temperature, in the rescaled specific heat 
is observed for < q < 1 (not for q > 1) whenever the energy spectrum includes a 
positive portion. Although no evidence for phase transitions has been obtained for 
q > 1, these should not be excluded without further studies. Indeed, the specific heat 
critical exponent a is herein shown to be positive for q < 1, and it is known to be zero 
for q — 1 (Mean- field Approach). Consequently, it could well be that it is negative 
for q > 1, thus producing a soft thermal dependence of the specific heat (as herein 
observed!). 

C. Distribution Function of a Paramagnetic System 

In [27], Biiyiikkihc and Demirhan have attempted to establish a similarity between the 
random walk problem and a paramagnetic system such a way as to the distance covered in 
the random walk here corresponds to the total magnetic moment of a paramagnetic system. 
Similar to the set of the points visited in the random walk, the statistical ensemble of the 
orientation of the magnetic moments in a paramagnetic system could be considered as a 
set of self similar elements, i.e. exhibiting a multifractal structure. When the Shannon 
entropy is maximized by taking the ordinary constraints, namely (i) normalization of the 
probability and (ii) the finiteness of the mean square total magnetic moment < M 2 >, 
the Levy distributions of the multifractal structures are not obtained. In this case, the 
maximum entropy formalism should be modified by taking another entropy definition instead 
of Shannon one and new constraints appropriate with the new entropy. Following the way 
used in [43], Buyiikkihg and Demirhan have used the integral form of the Tsallis entropy (in 
k units) 
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S, = (12) 

where p g (M)dM is the probability of the total magnetic moment of the ensemble to be 
in the interval M and M + dM ; with the appropriate constraints (i) normalization of the 
probability and (ii) the finiteness of (M 2 ) q which is the g-expectation value of M 2 in the sense 
of Eq.(3). By using the undetermined Lagrange multipliers method, it is straightforward 
to obtain the distribution function which requires to have the same asymptotic functional 
form with the Levy distribution, which behaves for large M as p(M) ~ M -1-7 (0 < 7 < 2). 
Hence, by equating the asymptotic solutions one ends up with 

In particular, the limiting cases provide an interval 1 < q < 3 as it has been given in [43]. 

D. Single-site Callen Identity 

In a recent work [29] , Sarmento has taken into account of single-site Callen identity from 
the GTT point of view, and accomplished to calculate the critical temperature of the Ising 
ferromagnet within mean field approximation. When the Hamiltonian in Eq. (9) is used, it is 
easy to obtain the standard single-site Callen identity (i.e. the standard expectation value 
of the spin variable at the lattice site) such that : 

(s^ = (tanh(/^)> (14) 

where Ei ~ J2j -hj s j- By approximating the thermal average of the hyperbolic tangent by 
the hyperbolic tangent of the thermal average, standard mean field approximation can be 
obtained : 



(sj) 1 ~ tanh 
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(15) 



The derivation of the generalized single-site Callen identity is started by seperating the Ising 
Hamiltonian into two parts 



H = Hi + H 



(16) 



where Hi includes all contributions associated with the site i and H' is the other part which 
does not depend on site i. After some algebra, within GTT, one can find 



l + f q 



(17) 



where 



l-P(l-q)(E i + H') " q/(1 - q) 



(18) 



l-/3(l- q )(-Ei + H') 

which is the generalized single-site Callen identity. After establishing this expression, he 
focuses his attention to the calculation of the critical temperature T c for a ^-coordinated 
spin-1/2 Ising ferromagnet with coupling constant J. In the T — > T c limit, the expression 
for the critical temperature is found to be 



kTc 
J 



= qz 



(19) 



which transforms to the well known standard result fcsT c /J = z in the q — > 1 limit. 



E. Mean-field Ising Model 

The authors of the present paper have recently developed a mean-field approximation 
to the Ising model within GTT [30,31]. In their study, they have used the generalized 
Bogolyubov inequality obtained in [25]. In the calculations, a trial Hamiltonian 

#o = -A]>> (20) 

i 

which corresponds to a system composed of N non-interacting single spins, is considered. 
From the partition function corresponding to a system with just one spin, the free energy of 
the N non-interacting spins system has been obtained by making use of approximating the 
partition function of the system with iV independent spins as the product of iV partition 
functions for systems with single spin. This "factorization scheme" is the same as the one 
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used in [22] to obtain the quantal distribution functions. By the help of this approximation, 
after some algebra, generalized mean-field magnetization and generalized mean-field free 
energy have been established respectively : 



l + PJz(l-q)((s) 



0)q 



PJz(l - q) ((s) 



9/(1-9) 



l+PJz(l-q)(( 8 ) ) 



9/(1-9) 



+ 



l-(3Jz(l-q)((s) ) 



9/(1-9) 



(21) 



[F mf ) c 



(Fo)Jl + (3(1 - q)±JzN((s) 2 ) Q ] + \JzN{(s)l\ 



l + P(l-q)JzN(( S y Q \ 



where 



(Fo) 



On - 



{[l + PJz(l-q)((8) \ 



1/(1-9) 



+ 



l-PJz(l-q)((s) ) q ] /{ 9) } -1 



(22) 



(23) 



0(1 -q) 

By investigating the graphical solutions of the magnetization and determining the minima 
of the free energy for various values of q, an interval has been established for the Tsallis q 
index. This result is found to be the same as the interval obtained for the paramagnetic 
free spin systems in [27], namely 1 < q < 3. The other remarkable derivation of this work is 
that the critical temperature of the generalized mean-field Ising model has been evaluated 
such that : 



kJc 
J 



= qz 



(24) 



which coincides with the results calculated in [29,39]. 



A few magnetic systems, which have been investigated within GTT, are also present, 
however we are contented to indicate what the subjects of these papers are about : 
1) For the spin-1/2 Ising ferromagnet, a transmissivity variable which extends that defined 
for thermal magnetic systems has been proposed. By using this generalized transmissivity 
as well as duality arguments, the (/-dependence of the critical temperature corresponding to 
the square lattice has been calculated [38]. 
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2) In [39], the authors have calculated the phase diagram and the correlation length critical 
exponent v for the Ising ferromagnet in a self-dual hierarchical lattice which mimics the 
square lattice. 

3) In order to illustrate the existance of a well behaved thermodynamical limit for the 
canonical ensemble, the ideal paramagnet has been discussed numerically. In addition to this, 
generalized Schottsky anomaly and generalized Curie law are calculated and data collapse 
is exhibited for iV >> 1 [40]. 
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